Roughness with a finite correlation length in the Microtrap 
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We analyze the effects of roughness in the magnitude of the magnetic field produced by a current carrying 
microwire, which is caused by geometric fluctuation of the edge of wire. The relation between the fluctuation 
of the trapping potential and the height that atom trap lies above the wire is consistent with the experimental 
data very well, when the colored noise with a finite correlation length is considered. On this basis, we gen- 
erate the random potential and get the density distribution of the BEC atoms by solving the Gross-Pitaevskii 
equation, which coincides well with the experimental image, especially in the number of fragmentations. The 
results help us further understand the nature of the fluctuation and predict the possible application in the precise 
measurement. 
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I. INTRODUCTION 

Trapping and manipulating ultracold atoms and Bose- 
Einstein condensations (BEC) in magnetic potential produced 
by micro atom chip attracted more attention recently because 
it is easy to perform coherent manipulation, transport and in- 
terferometry on a chip |l|,|2|], where the condensate is very 
close to the surface of wire. It is possible to bring miniatur- 
ization into the application of BEC [3]. The trapping frequen- 
cies in the radial direction are normally several kHz, and a 
few Hz in the axial direction. Hence the aspect ratio is so high 
that the size of BEC is several mm in the axial direction, the 
spatial effects on the condensate will be obvious. However, 
for small-scale magnetic field of atom traps produced by mi- 
cro fabricated current-carrying wires, an unexpected problem 
in these chips was their small fluctuation in magnetic field, 
which introduces a random potential along the trap. Such a 
fluctuation in potential will cause fragmentation of the atom 
clouds J3, 0] ■ The effects of the random potential have been 
studied theoretically [5, 6], where the randomness in magnetic 
field is due to geometric fluctuation of the wire surface and 
the strength of the disorder field shows a d~ 5 ' 2 relation with 
the height d that atoms lie above the wire, when assuming the 
correlation of the wire edge is the white-noise correlation. It 
is not coincidental to the experimentally observed d~ 22 rela- 
tion Q]. 

The deviation in the theoretical calculations JH from the 
experiment is caused by the white-noise approximation made, 
which does not meet the actual situation that the system's fluc- 
tuation cannot be the ideal white noise lHU. In this paper we 
elucidate the nature of the random potential, indicating the 
relation between the strength of the potential and the height 
of the BEC atoms. Different from the white-noise assump- 
tion, we consider a more general correlation form — the col- 
ored noise correlation R(r) = cr 2 e~' r ''^, with £ the correlation 
length of the system. Here we borrow the notion 'colored 
noise' from the theory for time-domain noise analysis to deal 
with the case of spatial roughness ifloll . Our theory shows that 




FIG. 1: (Color online) A typical atomchip for BEC. A constant cur- 
rent I flows in the z-direction of the wire with H high and Wo wide. 
The BEC atoms lie above the wire at the height of d. A bias magnetic 
field Bbias in the longitudinal direction produces a local minimum. 
/s(z) is the deviation of the edge from its ideal position. 



'Electronic address: xjzhou@pku.edu.cn 



the power index a in the disorder potential u s = C ■ d~ a is rele- 
vant to the system's intrinsic correlation length This means 
by measuring the dependence factor a, the correlation length 
of the current-carrying wire can be known. We point out in 
experiment 0], a = 2.2 because £ = 75/im. Consequently, 
the correlation function for the random potential is calculated. 

With the random potential we solve the Gross-Pitaevskii 
equation, depict the density profile of the captured BEC atoms 
and compare it with that of the white-noise case. The de- 
sired random potential is generated using a selected aperture 
function derived from the correlation function. The number 
of fragmentations has a better agreement with the experiment 
than previous calculations. This work is meaningful in the 
simulation of transport process. It as well predicts a criterion 
forjudging the quality of the wire used for cold-atom manip- 
ulation, which might be highly applicable in precise measure- 
ment. 
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H. THE STRENGHTH OF THE DISORDER POTENTIAL 

A typical setup of a magnetic micro trap is shown in Fig. 1. 
The current Iq flowing in a micro-fabricated copper conduc- 
tor in the z-direction and a bias field Bbias m the y-direction 
form a magnetic trap for atoms. The distance d from the 
surface of the current-carrying wire to the center of the trap- 
ping potential is controlled by the bias field and the current, 
d = idoIo/27iB bias . 

The random potential is brought about by the roughness in 
the magnetic trap for the BEC atoms lying above the wire sur- 
face about lOOyura ~ 200jura. What we are interested is the 
z-direction of the random magnetic field B z , which connects 
with the y-component of the current. In experiments 171 1 1 ill , 
the wire normally is a thin electroplated gold layer of a few yon 
thick. The surface of the wire can be made very flat techni- 
cally with the electroplating method. So the main contribution 
of the random field arises from the side of the wire. 

The distortion of the side edge of the wire curves the current 
flow. The potential roughness in the magnetic microtrap is in- 
duced by this abnormal current. To compute the fluctuation 
without loss of generality, we make the following assump- 
tions. First the current Iq remains constant on the edge of the 
wire, the amplitude Jq of the current density j is then con- 
stant too. Second the resistivity keeps unchanged along the 
current-carrying wire. Third the deviation of the edge from 
its ideal position (both left and right) /l/r is trivial compared 
with the width of the wire Wo, /l/r/ ^ Wo. Considering sym- 
metry, only f(z) = \{fi{z) + /r(z)) contributes notably to the 
fluctuation, and f~ = \(Jl — /«) is neglected. 

If we ignore the change in the module of the current den- 
sity vector, then we consider the effects due to the alteration 
of the direction of j. The z-component of B, i.e. B z , is directly 
relevant to the y-component of current density. In yz plane 
the current density satisfies the charge conservation condition 
^ + % = 0. Under such assumptions, we introduce an aux- 
iliary scalar potential G, so that j y = /ofp an d jz - 
The function G satisfies the Laplace equation in the interior of 
the wire V 2 G = 0. Using separation of variables, we get the 
Fourier component of f(z) which brings a transverse current 
density |5||, 



d r 



e ikz dk 

jy Kj 01T I . . .,„■ ■cosh(ky)smh(kW /2)F f (1) 
sinh(A;Wo) 



where Ff is the inverse Fourier transform of f(z), Ff - 
j f dze' kz f(z), and k is the wave vector in the frequency do- 
main. 

Therefore the disorder magnetic field can be derived from 
the Bio-Savart law: 
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The random potential U(z) is given by U(z) = -f*B ' B z . So 
the strength of U induced by the roughness in the magnetic 



microtrap can be measured by the value at the zero point of 
the potential's auto-correlation function g(r) with r — z - z', 



g (r) = (6U(z)6U(z')) 



(3) 



(...) denotes the ensemble average. And the Fourier transform 
of g(r), i.e. the power spectral density of the random potential, 
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As the wire is eletroplated with a thin layer of gold fTlll . it is 
acceptable to assume Hq <k d and ignore the thickness of the 
line. We substitute the dimensionless variable B — kd for k , 
and obtain 



S(B)=Ad- A -ffPiB) 



(5) 
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where A is a constant, F(k) = J d(z - z')c" ftfe " z ' ) </(z)/(z')>, 
K n (x) the Bessel function of the second kind, and 
y(n, x\, xo) is the incomplete gamma function j{n,x\,xi) — 
£ 1 dxx"~ l e~ x . We then compute the inverse Fourier transform 

and get the disorder potential u 2 s . 



g(z) = T- 1 [S(J3)] 
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When the correlation of the wire edge fluctuation is white- 
noise form, (f{z)f(z')) is a 6 function, and F(k) in Ecl(8) is a 
constant, we have « 2 ~ d~ 5 / 2 , which is the same in ||5|]. The 
white noise condition when the edge of the wire totally not 
correlated is the simplified case, it limits the supremum of the 
system fluctuation. Below we discuss the colored-noise model 
with a finite correlation length. 



III. DEPENDENCE BETWEEN POTENTIAL AND 
HEIGHT FOR A FINITE CORRELATION LENGTH 

The fluctuation f(z) is described as a colored noise model. 
The function f(z) is a Gaussian random variable with zero 
mean (f(z)) = 0, and variance iHHl] 



</W(z')> = cr 2 e 



-\m 



(9) 



The colored noise given by Eq.© is parameterized by three 
key factors, the distance along the wire r - z - z', the char- 
acteristic correlation length ^ and the noise strength cr. When 
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f — » 0, it degenerates into the white-noise correlation form 
(J(z)f£z?) = S(z-z'), which is the ideal situation of the rough- 
ness IS la]. When £ — > oo, the surface of the wire becomes 
totally correlated, implying a perfect smooth surface, which is 
the ideal case of the wire. The power spectrum of Eq.fO for 
the wire edge is 
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FIG. 2: (Color online) The power spectral function of potential S(J3) 
for the different correlation length of disorder potential. The solid 
line is for white noise with £ = 0/im. Others are for the case of 
colored noise. The dashed line with £ = 30/im, the dotted line £ = 
75fim, and the dash dotted line £ = 200/jm. 

When it comes to our model that the wire edge has an intrin- 
sic randomness, the power spectral density of which described 
in Eq.dTOb is moderate in comparison. The power spectral 
density of potential S(J3) for different correlation lengths is 
shown in Fig. 2. In the white-noise correlation £ = Qfj.ni, the 
curve peaks at f3 = 1.3. If we choose the typical correlation 
length £ = 15 urn 0], the curve peaks at/? = 1.1. From Fig. 2, 
we know the power spectral density peaks at a smaller value 
/? with the increasing ^. This means wave vector k decreases 
and the periodicity weakens for fluctuation at a certain d. The 
height of the peak goes up as the correlation length increases, 
which shows qualitatively the wire edge becomes smooth so 
that the frequency region is smaller for disorder potential. 

In order to clarify the effect of the wire width Wo, we 
plot the relation between the random potential u s and the ra- 
tio d/Wo in logarithmic coordinates in Fig. 3. To express it 
clearly, we can divide the figure into three regions. For dis- 
tance d varying from to Wo, the BEC atoms are very close 
to the top surface of the wire. The relation between the ran- 
dom potential and the ratio is a linear dependence, which is 
the fat-line approximation. When Wo < d < 3 Wo, the situa- 
tion becomes complicated, the rim effect of the wire as well 
as the orbiting magnetic field combine to present a field that 
cannot be easily calculated theoretically. For distance d rang- 
ing from 3Wo to infinity, this is also a linear fit corresponding 
to the thin-line condition. The BEC atoms are far from the 
top, so the width of the wire can be neglected and the wire as 
a whole dominates the field. In common experimental con- 
ditions SdU], tne atoms are placed high above the surface 
of the wire with Wo neglected. So below, during our discus- 
sion on the a ~ £ relation, the parameter d we choose mainly 



ranges from 100/im to 200jum. 
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FIG. 3: (Color online) The scaled power spectral density of the ran- 
dom potential and the dimensionless factor d/ Wo in logarithmic co- 
ordinates. In the area of < d/Wo < 1 and d/Wo > 3, the curve has 
a good linear fit. 

The dependence between the strength of the random poten- 
tial and d can be assumed as 

u s = C(d)-d-° (11) 

Based on the above analysis, the distance that the atoms lie 
above the wire d have a linear fit in the logarithmic coordinates 
for the thin-line approximation when d » Wo , so C(d) is a 
constant describing the amplitude of the random fields, a the 
power index characterizing the speed the fields decays with d. 
The strength of the disorder potential drops with the growth 
of height d. 

Using Eq.® and Eq.(fTTT). we calculate u s for different £ 
and fit them against d to get a. The relationship between a 
and £ is plotted in Fig. 4. We can see the potential roughness 
decays with the height d in direct relation to a. As the cor- 
relation length £ increases, the index a decreases drastically. 
When £ is zero, a — > 2.5, which is coincidental to the white- 
noise approximation. Provided that the deviation between the 
data from the experiment and theoretical calculation is mainly 
caused by the uncertainty of then in our calculations we 
elucidate the correlation length of the system, i.e., the fluc- 
tuation of the edge of the wire should be a finite number on 
a large scale. For instance, the correlation length in exper- 
iments [HI] should be 75/im pointed out by our figure as 
well as further calculations. The power index a is determined 
as soon as £ is determined. 

The strength factor <x 2 in the colored-noise correlation func- 
tion Eq.(9) only contributes to the constant before the integral 
in Eq.(7) and is totally not relevant to the shape of the random 
potential, so it has no effect on the power index a. However 
we can see below its influence on the density profile is similar 
to the alteration of the distance d. 



IV. THE FORM OF ROUGH POTENTIAL AND 
FLUCUATIONS OF THE DISTRIBUTION 

We then discuss the disturbance of the distribution of BEC 
atoms caused by the random magnetic fields. The rough po- 
tential in the trap can be obtained in experiments by measuring 
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FIG. 4: (Color online) Relationship between the power index a and 
the correlation length f in the case of thin-line approximation (d » 
Wo). The cross point is the typical value f = 75fim and a = 2.2. 



the density profile of cold atom clouds via absorption imaging 
method. As the BEC atoms expand mainly in the z-direction, 
and are very tight in the other two dimensions, the distribu- 
tion can be described by the one-dimension Gross-Pitaevskii 
Equation(GPE) fH 



ih 



2m 



v 2v p + v(zy¥ + nuq m 2 v 



(12) 



*F satisfies the unitary normalization condition, N is the num- 
ber of BEC atoms and Uq = 4nfi 2 a/m is the intensity of mean 
field interaction. V(z) includes both harmonic potential and 
the fluctuated potential as follows lfl3ll 



1 



V{z) = -mojh 2 + V d (z). 



(13) 



Vd(z) is the random potential induced by the fluctuation, and 
(x> z is the frequency in the z-direction. Here we use a method 
similar to IU4I1 to generate a random potential with desired 
correlation function. 



V d (z) = T- 1 [W(k)r[g(z)]], 



(14) 



where W(k) is a an appropriate aperture function which can 
filter out high and low part of the signal. With Vd(z) we use 
backward euler finite difference (BEFD) to compute the dy- 
namic evolution of the GPE and finally get the wave function 
*F and the density profile |*F| 2 fll . 

Due to the random potential induced by the fluctuation of 
the magnetic fields, a series of fragmentations in the conden- 
sate density profiles come into being. For different bias field 
values, the experimental sequence used to produce ultra-cold 
atoms in this trap were carried out, corresponding to different 
distances from the magnetic trap to the top surface of the wire. 
The strength of the random potential is mainly affected by the 
value cr 2 in Eq.(9). When the strength is relatively large, the 
random potential is comparable to the harmonic potential. The 
fragments present a high peak, and the density distribution 
rises and falls sharply, as shown in Fig. 5. 

Fig. 6 shows effect of the disorder strength on the den- 
sity profile. We focus on the area around the main peak and 
take out the 1000;um area. The relative roughness strength 
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FIG. 5: (Color online) Density profile of the BEC atoms and the 
random potential at different heights on the case of large random 
potential. The left column indicates the random potential at 50/jm, 
100/jm, 150/jm and 200/im respectively from top to bottom. The right 
column is for the normalized distribution \T\ 2 at the height of 50fim, 
100;um, 150yum and 200/im from top to bottom. Figures for I*?! 2 range 
the same from to 0.2. f is fixed at ISjim. [7J] 



in Fig. 6 (a-d) is 1:2:6:10. From the figure we can see the 
fragmentation gradually becomes obvious as the strength in- 
creases. At the appropriate strength for the random potential, 
as shown in Fig. 6(c), the density profile basically coincides 
with the experiment j 7] . The amplitude for the fragments are 
different with one highest and others decreasing accordingly, 
and the distance between two peaks in the lOOOjura area is 
about 250fim. The characteristics of the experiment J3l are 
reflected by our calculations. 
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FIG. 6: (Color online) Effect of the disorder strength on the density 
profile with £ = 75//m. The relative strength in (a), (b), (c) and (d) is 
1:2:6:10. 

First we discuss the effects of height d at a finite correla- 
tion length ^. We can see the dependence between the ran- 
dom potential and the distance in the left column of Fig. 5. It 
indicates the random potential at 50yum, lOOji/m, 150/jm and 
2Q0fim. The right column is for the density distribution. Us- 
ing the equation u 2 = CdT 22 given above, we can assess 
the strength of the potential quantitatively. As the distance 
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increases, the roughness of the random potential weakens. 
Meanwhile the fragmentation phenomenon in the density dis- 
tribution becomes less obvious. 
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FIG. 7: (Color online) The normalized density profile of the BEC 
atoms for different correlation length. The left column for the height 
d = 50pm. The right column for d = 100pm. All the figures have 
the same range for IT*! 2 from to 0.2. (a) £ = 0/jm (b) £ = 75/im (c) 
f = 200pm. 

Then we show the influence of different correlation length 
at a fixed height d. Fig. 7 demonstrates how the correlation 
length of the system affects the distribution of BEC atoms. 
The left column is for different correlation lengths 0pm, 75pm, 
200pm at the height d = 50pm, and the right column is for all 
the same conditions except for d = 100pm. When the sys- 
tem has an intrinsic correlation length, the fragmentation of 
the density profile is apparently small than that of the white 
noise. In two extreme cases, as £ — > oo, the distribution be- 
comes smooth because it is the ideal situation that the current- 
carrying wire has no fluctuation, while £ — > is the sim- 
ple white-noise case, and the fluctuation is at its maximum. 
We are most interested in Fig. 6(c) when d = 100pm and 
£ = 75pm. This set of data is measured in the experiment 10]. 
The number of the fragmentations, which is about 4 along the 
z-axis over 1000pm, is in good agreement with the image de- 
rived from the experiment, better than that of the white-noise 
approximation (about 6) which is different from the experi- 
mental data. 



V. DISCUSSION AND CONCLUSION 

According to our analysis, different correlation length will 
lead to different a. As in experiment fiEl, a = 2.2, the plau- 



sible correlation length through our calculations, should be 
75pm, and the density profile we thus get is coincidental with 
the experiment. The spectral density of the distribution Hill 
shows that £ has a very small fluctuation 100«m, so the fea- 
tures of a system could be considered as a finite correlation 
length. We also take the width of the wire into account, point- 
ing out exactly how it affects the random potential. Though 
the real trap is not a thin-wire, our model indicates when the 
height d is relatively large, it is rational to ignore the width of 
the wire. We then analyze how the strength of the random po- 
tential affects the density distribution. The density profile for 
the appropriate disorder strenth in our theoretical picture is in 
better accordance than the previous white-noise approxima- 
tion to some extent, which proves our conclusions to be more 
reasonable. That is to say, by counting the number of frag- 
mentations of the atoms, we can see qualitatively the correla- 
tion length of the wire. By furthermore measuring the power 
spectral density of the disorder potential, we can get the cor- 
relation length quantitatively. 

In summary, we bring forth a more general model to cal- 
culate the random potential for the BEC atoms in a current- 
carrying wire. We show that the fluctuation in the potential 
is caused by the distortion of the wire edge with a fixed cor- 
relation length. The correlation length of the specific wire 
used in experiment [0] is £ = 75pm and the power index in 
u s = C ■ d~ a is 2.2. This might be applicable in the precise 
measurement because of the possible criterion provided for 
the necessary wire quality to be used for cold atom manip- 
ulation. Besides, by measuring the spectral density of cold 
atom distribution and analyzing the deviation in the data we 
get from theoretical calculations and experiment, we can get 
to know the nature of the noise, and apply it in the simulation 
of the transport process. 
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